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Abstract. We model the acceleration of electrons at a 
shock front in a relativistic blazar jet and compute the ra- 
diation they emit in a post-shock region which contains a 
homogeneous magnetic field. The full space, time and mo- 
mentum dependence of the electron distribution is used 
in this calculation. It is shown that the 'homogeneous' 
synchrotron model is recovered, provided the downstream 
speed of the plasma away from the shock front is nonrel- 
ativistic, and provided that the light travel times across 
the face of the shock front is unimportant. By varying 
the rate at which particles are picked up by the accel- 
eration process, we calculate the time-dependence of the 
spectra. Since the magnetic field strength is assumed con- 
stant within the emission region, each frequency band can 
be identified with electrons of a particular energy. We 
find that for a band in which the electrons are acceler- 
ated rapidly compared to the rate at which they cool, the 
spectra typically harden during phases of rising flux, and 
soften during phases of falling flux, as has been observed 
in the objects PKS 2155-304 and Mkn 421. However, in 
a frequency band in which the timescales are compara- 
ble, the reverse behaviour is to be expected. We discuss 
the extent to which observations of both the stationary 
spectrum and the spectral variability of the synchrotron 
component of blazar emission can be used to constrain the 
model. 
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1. Introduction 

The nonthermal, highly variable X-ray emission of blazars 
is usually interpreted as the synchrotron emission of rela- 
tivistic electrons accelerated in a jet which itself moves at 



relativistic speed towards the observer. In this scenario, 
the higher energy gamma-rays detected from several such 
sources arise from the inverse Compton scattering of sof t 
photons by these electrons (for a review, see Sikora 1994 ). 
Two different types of model of the emitting region can 
be found in the literature: the 'inhomogeneous' and the 
'homogeneous' models. In the inhomogeneous model, the 
emitting part of the jet is assumed to be cylindrically sym- 
metric with a cross-section which varies with distance from 
the central object. The magnetic field is also taken to vary 
with distance from the central object. It is usually as- 
sumed that relativistic electrons are 'injected' at a shock 
front which moves through the emission region, starting 
close to the central object. These electrons also move out- 
wards with the jet material, and undergo energy losses due 
both to synchrotron radiation and to the adiabatic expan- 
sion of the jet (Marscher & Gear 1996| , Maraschi, Ghis- 
ellini & Celotti |1992| , Marscher & Travis |1996[ ). The re- 
sulting synchrotron spectrum is found by integrating over 
the emission region, within which both the magnetic field 
and the particle distribution is inhomogeneous. 

In the homogeneous model, on the other hand, both 
the magnetic field and the particle distribution function 
are assumed homogeneous throughout the emission region 
(Inoue & Takahara ^996|, Chiaberge & Ghisellini ^997 
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Mastichiadis & Kirk 1997). The relativistic electrons are 
injected with a specified distribution and are assumed to 
escape on a timescale i csc . After escape, a particle no 
longer radiates. Although at first sight somewhat arbi- 
trary, there is a good observational reason to believe that 
escape, or, equivalently, sudden energy loss by adiabatic 
expansion, is important: the spectral index of the radio 
emission of blazars is typically hard (a > —0.5) and can- 
not be produced by particles which have been allowed to 
cool completely by synchrotron emission. If we adopt a 
strictly homogeneous model, electrons must escape from 
the emission region and thus be prevented from cooling 
completely. Such a picture departs from the standard ex- 
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planation, in which the flat spectrum is thought to be the 
result of a variation of the self-absorption frequency within 
the source in an inhomogeneous model (the 'Cosmic Con- 
spiracy' Marscher 1980). Depending on the magnetic field 



strength and Doppler boosting factor, self-absorption may 
also become important in a homogeneous model. Addi- 
tional components emitting at low frequency would then 
be required. 

The homogeneous model with escape corresponds ap- 
proximately to the plausible physical situation in which 
particles are accelerated at a shock front, provided there 
is a region of relatively high magnetic field just behind 
the shock. In this case, radiation from this region may 
be expected to dominate the observed emission. This re- 
quires that on leaving the region (on a timescale i esc ), 
the particles encounter a magnetic field that is so weak, 
that no significant further contribution to the emitted 
spectrum arises, despite the fact that particles accumu- 
late there over the entire life of the source. However, a 
problem arises in the homogeneous model concerning time 
variability - only if the physical dimensions of the source 
are such that the light crossing time is short compared 
with the synchrotron cooling time is it reasonable to as- 
sume homogeneity. If this condition is not fulfilled, the 
observed variability is dominated by the shape and ori- 
entation of the source, rather than the intrinsic cooling 
and acceleration timescales. Until recently, models of the 
rapid variability of blazars have either fulfilled this condi- 
tion or have made specific assu mptio ns about the source 
geom etry (Mastichiadis & Kirk 1997 , Chiaberge & Ghis- 
ellini |l997|) . 

In this paper we present a model in which particles 
are accelerated at a shock front and cool by synchrotron 
radiation in the homogeneous magnetic field behind it. 
The plasma downstream of the shock front moves rela- 
tivistically towards the observer; the shock front is non- 
relativistic when seen from the rest frame of the down- 
stream plasma. The kinetic equations are solved for the 
time, space and energy dependences of the particle distri- 
bution function, and the resulting synchrotron emission is 
calculated. 'Escape' is accounted for by assuming that the 
magnetic field strength drops suddenly at a finite distance 
behind the shock, so that the radiation from the adiabati- 
cally cooled electrons in the weaker field can be neglected. 
Thus, this model is homogeneous in the sense that the 
magnetic field does not vary through the emission region, 
but contains an inhomogeneous electron distribution. The 
variability is computed by assuming the observer lies in 
the direction of the normal to the shock front. This pre- 
ferred orientation does not affect the observed variability 
provided the light travel time across the face of the source 
(i.e., over the surface of the shock) is short compared to 
the synchrotron cooling time as measured in the rest frame 
of the plasma. At the highest electron energies considered 
in the application to X-ray blazars, this condition may be 
violated, in which case the computed variability will be 



smoothed out over the longer timescale. Whether or not 
this is expected to occur could be decided by comparing 
the predicted emission by inverse compton scattering with 
gamma-ray observations. In principle, observations of the 
frequency dependence of the variation timescale could also 
resolve the question. 

We consider only the synchrotron radiation of the ac- 
celerated particles, leaving the more involved computation 
of the inverse Compton emission to future work. This is 
sufficient for comparison with the observed radio to X-ray 
emission of blazars, provided the energy losses of the elec- 
trons are not dominated by inverse Compton scattering, 
which is usually the c ase, at least for BL Lac-type objects 
(Comastri et al. |1997| ). 



2. The electron distribution function 

Consider, then a shock front propagating along a cylin- 
drical jet of constant cross-section. Electrons are acceler- 
ated at the shock front, and subsequently drift away from 



it in the downstream flow. Following Ball & Kirk (1992), 
we treat two spatial zones: one around the shock front, 
in which particles are continuously accelerated, and one 
downstream of it, in which particles emit most of the ra- 
diation. In each zone, it is assumed that pitch-angle scat- 
tering maintains an almost isotropic particle distribution. 
The number iV(7)d7 of particles with Lorentz factor be- 
tween 7 and 7 + d7 in the acceleration zone around the 
shock is governed by the equation 



ON 
~dt 
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d~f 
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+ = 0*(7"7d) (1) 
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(Kirk et al. |199J), where 

4 cr T f B 2 
3 m c c \ 2 no , 

with or = 6.65 • 10 _29 m 2 the Thomson cross section. The 
first term in brackets in Eq. (1) describes acceleration at 
the rate the second describes the rate of energy loss 
due to synchrotron radiation averaged over pitch-angle 
(because of the isotropy of the distribution) in a mag- 
netic field B (in Tesla). Particles are assumed to escape 
from this region at an energy independent rate t~ s ]., and to 
be picked up (injected) into the acceleration process with 
Lorentz factor 70 at a rate Q particles per second. Note 
that the concept of this 'acceleration zone' differs from the 
emission region in the usual homogeneous models in two 
important respects: a) particles are injected at low energy 
and continuously accelerated and b) very little radiation 
is emitted by a particle whilst in the acceleration zone (see 
the discussion below). 

Assuming a constant injection rate Qo after switch-on 
at time t = 0, and setting N(j) = for t < 0, Eq. (1) has 
the solution 

(tacc-*esc)/te 

N(j,t) = a- 
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7n 



J.G. Kirk, F.M. Rieger & A. Mastichiadis: Particle acceleration and synchrotron emission in blazar jets 



3 



x0(7-7o)6(7i(t) -7), 



(3) 



for 70 < 7 < 7i(t), and N(j, t) = 0, otherwise. Here Q(x) 
is the Heaviside step function, 



a — Qat acc ^ acc /*° 



70 



7n 



-tacc/tc 



and the upper bound 71 (t) is determined by 



7i W = 



(S + 
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(4) 



(5) 



with 7 max = (/3 s i aC c) 1 . The quantities t acc and t esc have 
been assumed independent of energy in this solution, the 
more general solution allowing for energy dependence is 
given in the Appendix. 




Fig. 1. The electron distribution integrated over the source, 
as given by Eq. (10). The three curves correspond to the times 
5t a cc (solid line), 3ut acc (dotted line) and 500t acc , (dashed line). 
At larger times the distribution does not evolve appreciably. 
For this plot, t csc — 2t acc (i.e., s = 3.5) and tb = 100i acc . 
At t = 5t a cc, no particles have had time to cool, since the 
cooling time at the maximum Lorentz factor of 7 ~ 0.0257 ma x 
is approximately 40t acc . At times t > 500£ acc all particles with 
7 > 7maxtacc/ib cool, whereas those of lower 7 leave the source 
without significant loss of energy. 



In the model developed by Ball & Kirk (1992), ac- 
celerated particles escape into the downstream plasma, 
where they radiate. We can formulate the kinetic equa- 
tion obeyed by the density of particles in the radiation 
zone most compactly using a coordinate system at rest 
in the radiating plasma. The shock front then provides 
a moving source of electrons, which subsequently suffer 



energy losses, but are assumed not to be transported in 
space. The kinetic equation governing the differential den- 
sity dn(x, 7, t) of particles in the range dx, dj is then 



at c»7 i esc 



(6) 



where x s (t) is the position of the shock front at time t. For 
a shock which starts to accelerate (and therefore 'inject') 
particles at time t = and position x — and moves at 
constant speed u s , the solution of Eq. (6) for 7 > 70 is 



1 



- A * - 



1 



-, (*acc-*esc)/*e 



7 V u sJ 7max. 

e[ 7l (x/ Us )-(l/ 7 -/3 s t + / 3 s x/ Us )- 1 ] , (7) 



where 71 (t) is given by Eq. (5). To obtain the synchrotron 
emissivity as a function of position, time and frequency we 
convolve the density n with the synchrotron Green's func- 
tion P(v, 7). A convenient approximation to this function 
is given by Melrose fll98C| ): 



P{y, 7) = a s z exp(— z) 



(8) 



where a s = \^3e 2 fl sin0/(27rc) is a constant and z = 
47ri//(3f2sin(97 2 ), with £1 = eB/m the electron gyro fre- 
quency and 9 the angle between the magnetic field di- 
rection and the line of sight. At a point x — X (> u s t) 
on the symmetry axis of the source at time t the specific 
intensity of radiation in the x direction depends on the 
retarded time t = t — X/c and is given by 



(9) 



h{v,t) = / d7P(^, 7) / dxn(x,j,t + x/c) 



and the integrated particle density can readily be evalu- 
ated: 



dx n(x, 7, t + x/c) 



(1 - u s /c)j 



max 



0/t. 



7 



t x(l — u s /c) 



7 




.(10) 



The limits of the spatial integration are given by the re- 
tarded position of the shock front 



Xi 



(*) = 



ud 



1 - (u s /c) 



(11) 



and the retarded position 2:0(7, t) of the point furthest 
from the shock front at which particles have Lorentz factor 
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7 at time t + xq/c. This is given either by the solution 
2 coo i(7,i) of the transcendental equation 



t -f~ X C ool/ C ^-cool 



1 



To 



1 ) exp [-x cool /(u s t 

acc )\ 



(12) 



or by the assumed maximum spatial extent of the emission 
region (i.e., the point at which the magnetic field declines 
substantially). Denoting this distance by L, we have 



x ("f,t) = Max[a; CO oi(7,t),:Ei(i) - L] 



(13) 



This is most conveniently expressed in terms of the (re- 
tarded) time t b for the plasma to traverse the emitting 
region, as measured in the plasma rest frame: 



(1 - u s /c)L 



(14) 



Equation (10) gives the integrated particle density for 
times greater than the 'switch-on' time: 



t > t c 



t a 



(7max/7o) - 1 
(7max/7) - 1 



(15) 



before which it vanishes. The resulting electron spectrum, 
integrated over the source is depicted in Fig. [|. A char- 
acteristic break in the spectral slope appears at a par- 
ticular Lorentz factor 7 br (i) which, at any given time, 
separates those electrons which cool within the source 
(7 > 7break) from those which do not cool within the 
source (7 < 7 brca k)- At large times, all electrons with 
7 < 7brcak(oo) = j max /(l + t h /t acc ) leave the source be- 
fore cooling, and the integrated electron density becomes 
time-independent. 

A question which remains open in this approach is 
the synchrotron radiation emitted by a particle whilst in 
the acceleration region. If the magnetic field ahead of the 
shock were the same as that behind the shock, the to- 
tal emission could easily be computed using Eq. (3). For 
t > ton, we would have an extra contribution to the flux: 

I 8 (v,t) = J d 7 P(*/,7) J dxN{j,i)5[x-u s (t + x/c)] 

2 



(1 - Vch^x +t ™>/*~ 



dryP(v, 7) 



'"/max 

7 



(16) 



and the total emission Ii{y,t) would be given by 
Jl(l/,t) = J (f,t)+J.(i/,t) (17) 

However, for oblique shocks, the magnetic field 
strength is expected to increase upon compression at the 



shock. Particles undergoing acceleration spend part of 
the time in the upstream and part in the downstream 
plasma, so that it is not clear how to evaluate their syn- 
chrotron emission, although Eq. (16) certainly gives an up- 
per limit. At oblique shocks, reflection at the front itself is 
thought to be more important than diffusion in the down- 
stream zone, (Kirk & Heavens 198S| ), so that accelerating 
particles spend all their time upstream. In this case, it 
seems reasonable to neglect the emission from the accel- 
eration zone completely, which is the approach adopted 
here. An improved treatment of this point demands a full 
time-dependent solution of the diffusion advection equa- 
tion, which requires considerable numerical effort (Fritz & 
Webb [L990D , 

All quantities calculated so far in this section refer to 
the frame in which the radiating plasma is at rest (the 
jet frame). For application to blazars, they must be trans- 
formed into the observer's frame. Assuming the observer 
to lie in the direction of motion of the plasma, and denot- 
ing the plasma bulk velocity by c(3 D , the relevant transfor- 
mations are 



X 
t 

t'obs 

Jobs(^obsi ^obs) 



r b (Z}-/3 b ci obs ) (18) 

r b (tobs - 0bD/c) (19) 

r b (i + A>«2r bi , (20) 

{v ohs / V fl{v,t) a &Tll{v,t) (21) 



where T b (^> 1) is the Lorentz factor of the approaching 
jet plasma and D the position of the detector with respect 
to the position of the shock at t = 0, as measured in the 
observer's reference frame. 

Several simple qualitative results follow from these ex- 
pressions. Close to the maximum emitted frequency, the 
timescale on which the intensity varies in the frame of the 
plasma is roughly the switch-on time ~ (1 — u s /c)t acc . Ac- 
cording to Eqs. (18) and (19), t = (1 + /3 b )r b (i obs - D/c) 
Thus, the observed timescale is shorter than the intrinsic 
iacc by a factor r b (l+/3 b )/(l — u s /c). Doppler boosting of 
the flux, is independent of the shock speed, and is given 
simply by Eq. (21). 

3. Application to the synchrotron spectra of 
blazars 

The simple model described above provides a remarkably 
good fit to the radio to X-ray spectra of several blazars. As 
an example, we show in Fig. || observations of the object 
Mkn 501. The gamma-ray emission of this object, which 
has been the subject of much recent interest (e.g., Brad- 



bury et al 1997), is not included in this figure, since it 



is not thought to arise as synchrotron radiation. For the 
X-ray emission we di splay the archival data selected by 
Catanese et al. ( 1997 ) and not the data taken during the 
TcV flare in April 1997. 

The stationary emission (t — ► 00) found from Eq. (9) 
by a single numerical quadrature over the synchrotron 
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Fig. 2. The radio - X-ray spectrum of the o bject Mkn 501 
(data taken from the collation of Catanese et al. 1997) together 



with the stationary synchrotron emission from a single homo- 
geneous source 

Green's function is also shown in Fig. ^. The form of 
the spectrum is very close to that given by Meisenheimer 
& Heavens ( 1987 ), who used an analytic solution to the 
stationary diffusion/advection equation, including syn- 
chrotron losses. Four free parameters are used to produce 
this fit: 

1. the low frequency spectral index a = —0.25, which 
corresponds to taking i acc = t csc /2 

2. the characteristic synchrotron frequency emitted by an 
electron of the maximum Lorentz factor as seen in the 
observers frame (taken to be 1.3 x 10 18 Hz) 

3. the spatial extent of the emitting region, which de- 
termines the position of the spectral break at roughly 
5 x 10 12 Hz, corresponding to = 700t acc 

4. the absolute flux level. 

Since we restrict our model to the synchrotron emis- 
sion of the accelerated particles, it is not possible indepen- 
dently to constrain quantities such as the Doppler boost- 
ing factor, or the magnetic field. These can, however, be 
found using a model for the gamma-ray emission, for ex- 
ample the synchrotron self-compton model (Mastichiadis 
& Kirk 1997). Similarly, the frequency below which syn- 



chrotron self-absorption modifies the optically thin spec- 
trum is not constrained within the current picture. Never- 
theless, our model of the synchrotron emission makes pre- 
dictions concerning the spectral variability in each of the 
three characteristic frequency ranges which can be identi- 
fied in Fig. ||. These ranges are generic features of any syn- 
chrotron model, so that the predicted variability can easily 
be applied to the synchrotron emission of other blazars. 



They are a) the low frequency region, where the particles 
have not had time to cool before leaving the source (this 
is the region with a = —0.25 in Fig. |[ below the break 
at 5 x 10 12 Hz) b) the region between the break and the 
maximum flux, where the particles have had time to cool, 
but where the cooling rate is always much slower than the 
acceleration rate and the spectrum is close to a — —0.75, 
and c) the region around and above the flux maximum at 
roughly 10 17 Hz, where the acceleration rate is comparable 
to the cooling rate. 

Variability or flaring behaviour can arise for a num- 
ber of reasons. When the shock front overruns a region 
in the jet in which the local plasma density is enhanced, 
the number of particles picked up and injected into the 
acceleration process might be expected to increase. In ad- 
dition, if the density change is associated with a change 
in the magnetic field strength, the acceleration timescale 
might also change, and, hence, the maximum frequency of 
the emitted synchrotron radiation. Considering the case 
in which the acceleration timescale remains constant, it is 
a simple matter to compute the emission, since Eq. (6) is 
linear. An increase of the injection rate by a factor 1 + r/f 
for a time tf is found by setting 



Q(t) — Qo for t < and t > t{ 
Q(t) = (l + r/ f )Qo for0<i <t { 

We then have 

t) — I\{u, oo) + 
m [h(y,i) - h{v,i- (1 - u s /c)t { )] 



(22) 
(23) 



(24) 



Using r/i — 1, tf = 10i acc and u s = c/10, we show 
the resulting emission at a frequency v = z/ max /100 in 
Fig. ||. In the case of Mkn 501, this corresponds to a fre- 
quency of about 10 16 Hz, which lies towards the high fre- 
quency part of region b), between the infra-red and X-ray 
regions, where the spectral index is close to a = —0.75. 
Also shown in this figure is the temporal behaviour of the 
spectral index, as determined from the ratio of the fluxes 
at 0.01^ max and 0.05f max , through the flare. When plotted 
against the flux at the lower frequency, the spectral index 
exhibits a characteristic loop- like pattern, which is tracked 
in the clockwise sense by the system. This type of pattern 
is well-known and has been observed at different wave- 
lengths in several sources e.g., OJ 2 87 (Gear et al. |1986|) , 
PKS 2155- 304 (S embay et al. [1993D and Mkn 421 (Taka- 
hashi et al. 1996 ). It arises whenever the spectral slope is 
controlled by synchrotron cooling, (or, in fact, any cooling 
process which is faster at higher energy) so that informa- 
tion about changes in the injection propagates from high 
to low energies (Tashiro et al. 1995 ). 

If the system is observed closer to the maximum fre- 
quency, where the cooling and acceleration times are 
equal, the picture changes. Here information about the oc- 
currence of a flare propagates from lower to higher energy, 
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Fig. 3. The intensity and spectral index during the flare de- 
scribed by Eq. (23), as a function of time at low frequency. 
The loop in the a vs. intensity plot is followed in the clockwise 
direction. 



as particles are gradually accelerated into the radiating 
window. Such behaviour is depicted in Fig. |[ where the 
same flare is shown at frequencies which are an order of 
magnitude higher than in Fig. ||. In the case of Mkn 501, 
the frequency range is close to 10 18 Hz. This time the loop 
is traced anticlockwise. Such behaviour, although not as 
common, has also occasionally been observe d, for e xample 
in the case of PKS 2155-304 (Sembay et al. |l993| ). 

Finally, in the region of the spectrum below the break, 
where the spectral index in the case of Mkn 501 is close to 
a = —0.25, the emission is determined by the finite size 
of the source region, as expressed by tb- Here, the flare 
shown in Figs. || and |] has only a very small effect on 
the observed flux, since the duration tf = 10t acc has been 
chosen to be much smaller than the time tb taken to fill 
up the emitting region with radiating particles. However, 
even in the case of a larger flare, or one of longer duration, 
no variation of the spectral index is to be expected through 
the flare at frequencies below the break, because the time 
taken to fill the effective emitting region is independent of 
frequency 

This effect of smaller changes in the spectral slope at 
lower frequencies is also evident from Fig. ^, where the 
spectrum is shown at times t = (i.e., the stationary 
emission as shown in Fig. ^|), t = 10£ acc , t = 20i acc , and 
t = 30t acc . Here it can be seen that the rise in emission is 
rapid at all frequencies, and the subsequent fall sets in as 
a wave which propagates downwards in frequency. 

The results presented in Figs. ||and|4|are computed for 
u s = 1/10, so that the effects of the finite light travel time 




10 20 30 40 50 
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Fig. 4. The intensity and spectral index in the same flare as 
in Fig. ^| but at high frequency. The loop in the a vs. intensity 
plot is followed in the anticlockwise direction. 
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Fig. 5. The spectrum at times t = (solid line), 10 (dot- 
ted), 20 (short-dashed) and 30t acc (long-dashed) during the 
flare shown in Figs. ^| and ^. Data points are taken from 
Catanese et al. ( |l997| ). 
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between the front and back of the source are negligible. 
In fact, in the limit u s — > 0, Eqs. (7) to (24) reduce to the 
equations for a homogeneous source, provided the limit is 
taken keeping finite, i.e., allowing the maximum spatial 
extent of the source u s tb to vanish. 

4. Discussion 

It is generally accepted that the emission mechanism re- 
sponsible for the radio to optical emission of blazars is 
the synchrotron process. In objects similar to Mkn 501, 
which is discussed in Sect. 3, the synchrotron emission ex- 
tends up to X-ray energies. It has recently emerged that 
a model in which the radiating relativistic electrons re- 
side in a region of uniform magnetic field can provide a 
reasonable fit to observations not only of the synchrotron 
component, but also, via the inverse compton scattering 
of the synchrotron photons, to gamma-ray observations in 
the GeV to TeV range. In this paper, we have introduced 
a detailed model of the emission region including the pro- 
cesses of acceleration and synchrotron cooling. This should 
ultimately help us to assess how accurately a homogeneous 
region can reproduce the observed synchrotron emission. 
In combination with a computation of the inverse comp- 
ton emission, the physical conditions and the geometry of 
the emission zone can be constrained. 

The model used to describe acceleration is similar to 



that used by Ball & Kirk (1992) to describe the nonther- 
mal radio emission of SN1987A. It involves a phenomeno- 
logical division of the particle population into two groups: 
particles which are undergoing both acceleration and cool- 
ing and particles which merely cool. As a result, distribu- 
tions are found which extend to a maximum value of the 
Lorentz factor, and then cut-off abruptly. The synchrotron 
emission from such a distribution fits well with that from 
the source we have discussed in Sect. 3 - Mkn 501. How- 
ever, the form of the spectrum, especially in the hard X- 
ray region, depends on the shape of the cut-off of the elec- 
tron spectrum. In our model, the sharpest possible turn- 
over is produced. Inhomogeneities within the source or 
refinements of the acceleration model result in a broader 
turn-over. This indicates that in the case of Mkn 501 a 
homogeneous model of the region in which synchrotron 
X-rays are emitted is a good approximation. It also con- 
firms that a detailed model of the cut-off - which would 
necessarily involve additional parameters - is unneces- 
sary. In other blazars e.g., Mkn 421 as well as in other 
synchrotro n emit ting objects e.g., SNR1006 (Mastichiadis 
& de Jager 1996 ), a broader cut-off is indicated by the 
observations. Our current computations suggest that if 
the physical situation is similar to that in Mkn 501, the 
broadening of the cut-off should be attributed to inhomo- 
geneities within the source rather than an intrinsic prop- 
erty of the acceleration mechanism. 

The variability predicted for a homogeneous source de- 
pends on the frequency of observation and the param- 



eters of the acceleration mechanism - for relatively low 
frequency radiation, a characteristic pattern is produced, 
as pointed out by several groups. We have shown that 
closer to the maximum emitted frequency, this pattern 
should change. However, the observed variability depends 
not only on the intrinsic time-dependence, but also on the 
smoothing caused by light travel time delays across the 
source. We have described in detail the situation when the 
observer is positioned exactly on the axis of the source so 
that the emitting plasma moves directly towards the ob- 
server. This assumption is not as restrictive as it might 
at first sight appear, since the emission from a source in 
relativistic motion with Lorentz factor T is significantly 
boosted when viewed from a direction which makes an 
angle of less than ~ 1/r with the velocity. Intrinsic time 
variations within the source then appear shorter by a fac- 
tor r (see Sect. 2). These are also smoothed out on the 
timescale i?/(cT), where R is a typical dimension of the 
source perpendicular to the line of sight, provided that the 
angle between the source velocity and line of sight is of the 
order of 1/r. It is a good approximation to neglect this 
smoothing if the intrinsic synchrotron cooling time, mea- 
sured in the rest frame of the emitting plasma, is longer 
than the light travel time R/c. The same restriction ap- 
plies to the size of the source along the direction to the ob- 
server. Here, however, the intrinsic variations always dom- 
inate if u s <C c. The general formulae given in Sect. 2 are 
valid also for relativistic u a , but in the examples predicted, 
we have restricted ourselves to u s = c/10. Whether or not 
intrinsic variability dominates over light travel time effects 
depends on the frequency of observation. For Mkn 501, for 
example, we have found that the maximum timescale over 
which particles cool is determined by the position of the 
spectral break and is roughly 700 times the synchrotron 
cooling time at the maximum emitted frequency. Thus, de- 
pending on R, there may exist a critical frequency above 
which variations are smoothed out by light travel time ef- 
fects, but below which the results of Sect. 3 are valid. A 
detailed model of the inverse Compton emission is needed 
to estimate this frequency. 
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and where we have written 

£acc(7) = 7*a,0 t csc (j) = jt cfi (A7) 

Note that in this case the power-law index s, defined by 
N (x j~ s+2 is related to the timescales by s — 2 + t a o/< c ,o, 
whereas in the case of energy independent acceleration and 
escape s — 3 -\- t acc /t esc . 



A. Energy dependent acceleration and escape 



In general, one might expect the acceleration time t acc and 
escape time t csc introduced in Eq. (1) to be functions of 
particle energy. In this case the solution Eq. (3) is modi- 
fied. Defining 



/(7) 



7 



t acc {l) 

which is positive in the range of interest, one can write 
N 



(Al) 



ON d 

-W + JT /(7 )N] nr — n 



Q e(t)6(>y - 7o ) • (A2) 



The solution of this equation subject to the boundary con- 
dition iV( 7 , 0) = is easily found using Laplace trans- 
forms: 



N( 7 ,t) = 



where 



Qo 

/(7) 



[t — t(7)] exp 



d 7 ' 



icsc(7')/(7') 



r( 7 ) 



dy 
7(7) 
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(A4) 



In addition to the straightforward case t acc = constant, 
tcsc — constant dealt with above, it is also interesting to 
consider the case in which both of these quantities are 
linearly proportional to 7. This would arise in modelling 
diffusive acceleration with a 'gyro-Bohm' spatial diffusion 
coefficient (e.g., Kirk et al. 1994). The solution is 



iV( 7 ,f) = 0[t-r( 7 )] 
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